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Particle acceleration timescales in relativistic shear flows 
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We review the acceleration of energetic particles in relativistic astrophysical jets characterized by a significant 
velocity shear. The possible formation of power-law momentum spectra is discussed and typical acceleration 
timescales are determined for a variety of different conditions such as parallel and azimuthal shear flows. Special 
attendance is given to the analysis of parallel shear flows with either a linear decreasing or a Gaussian-type 
velocity profile. It is shown that in the presence of a gradual shear flow and a particle mean free path scaling 
with the gyroradius, synchrotron radiation losses may no longer be able to stop the acceleration once it has 
started to work efficiently. Finally, the relevance of shear acceleration in small- and large-scale relativistic jets 
is addressed. 



1. Introduction 

Relativistic jet outflows are observed across a wide 
range of astrophysical scales, from galactic micro- 
quasars to GRBs and radio-loud AGNs 0, E E3|. 
The fact that velocity gradients are a generic feature 
of environments producing jets suggests that many 
(if not all) of these outflows may possess a signifi- 
cant internal velocity shear. Today, the phenomeno- 
logical evidence for internal jet stratification is indeed 
mounting and comprises evidence for universal struc- 
tured jets in GRBs [ill l23j . internal jet rotation in 
AGNs [2(], EI, and a multi-component jet structure 
consisting (at least) of a fast moving inner spine and 
a slower shear layer as indicated, for example, by the 
detailed analysis of the brightness and polarization 
systematics in kpc FR I jets 0, ^| or the inten- 
sityand polarization maps of pc-scale FR I + II jets 
[3l ll9l.l25| . It seems likely that in the presence of such 
a velocity shear efficient acceleration of particles to 
high energy may occur, provided particles are scat- 
tered across the shear by mag netic inhomogeneities 
carried within the flow [^j, |22| . The theory of shear 
acceleration deals with those cases where the impact 
of the systematic velocity components of the scatter- 
ing centers becomes dominant over their random mo- 
tion components known as the source of second-order 
Fermi acceleration effects. 

Substantial theoretical contributions in the field of 
shear acceleration have been given by several au- 
thors (cf. [2l| for a review): Berezhko & Krym- 
skii (1981), for example, showed that (non-relativistic 
gradual) shear acceleration can lead to power law par- 
ticle momentum distributions resembling those of the 
classical Fermi theory, whereas Earl, Jokipii & Mor- 
fi.ll (1988) derived the corresponding diffusive particle 
transport equation in the diffusion approximation in- 
cluding shear and inertial effects. The generalisation 
of the work by Earl et al. to the relativistic regime was 
achieved by Webb using a mixed-frame ap- 

proach and applied by Rieger & Mannheim (2002) to 
rotating and shearing jets in AGNs. Particle accelera- 
tion in non-gradual relativistic shear flows (i.e., in the 



presence of a relativistic velocity jump), on the other 
hand, has been analysed by Ostrowski 0, IT7L ITsf 
using Monte Carlo methods, showing that very flat 
momentum spectra may occur. 

2. Microscopic picture 

Shear acceleration draws on a simple mechanism 
which seems to be applicable to a wide range of astro- 
physical flows such as accretion flows and structured 
jets in GRBs and AGNs (cf. also |j)- According to 
the underlying physical picture particles gain energy 
by scattering off (small-scale) magnetic field irregu- 
larities with different local velocities due to being sys- 
tematically embedded in a collisionless velocity shear 
flow. The scattering process is assumed to occur in 
such a way that the particles are randomized in direc- 
tion, with their energies being conserved in the local 
commoving fluid frame. In the presence of a velocity 
shear, the momentum of a particle travelling across 
the shear will change with respect to the local scatter- 
ing frame so that for an isotropic particle distribution 
a net increase may occur Q . For illustration consider 
a continuous, non-relativistic shear flow with velocity 
field given by 

u = u z (x)e z . (1) 

Let v = (v x , v y , v z ) be the velocity vector, m the rela- 
tivistic mass and p\ the initial momentum (relative to 
local flow frame) of the particle. Within one scattering 
time t (initially assumed to be independent of momen- 
tum) the particle travels a distance Sx = v x r across 
the shear, so that in the presence of a gradual shear 
the flow velocity will have changed by Su — duel, with 
Su = (du z /dx) Sx. Hence the particle's momentum 
relative to the flow becomes p2 = Pi + m Su (Galilean 
transformation), i.e. 

p\ = p\ + 2m Su p 1:Z + m 2 (Su) 2 . (2) 

As the next scattering event preserves the magni- 
tude of the particle momentum (in the local scatter- 
ing frame) the particle magnitude will have this value 
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in the local flow frame. Using spherical coordinates 
and averaging over solid angle assuming an almost 
isotropic particle distribution, the resultant average 
rate of momentum change and the average rate of mo- 
mentum dispersion scale with the square of the flow 
velocity gradient jal23|, i.e., 
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~AT 
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du z 
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dx 
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Generalizing these results to a momentum-dependent 
scattering time obeying a power law of the form 
r oc p°, we may write down a simple Fokker-Planck 
transport equation for the phase space particle dis- 
tribution function f(p), assuming a mono-energetic 
injection of particles with momentum pq. Solving for 
the steady-state with a > one obtains (see [2^) 



/(p)cxp~( 3+Q > H(p-p ). 



(5) 



where H(p) is the Heaviside step function. Hence 
for a mean scattering time scaling with the gyro- 
radius (Bohm case), i.e. r oc p a , a = 1, one obtains 
f(p) oc p~ A , i.e., a power law particle number density 
n(p) oc p~ 2 which translates into a synchrotron emis- 
sivity j v oc v~ x l 2 . For a Kolmogorov-type (a — 1/3) 
or Kraichnan-type (a = 1/2) scaling, on the other 
hand, much flatter spectra may be obtained, i.e., 
the synchrotron emissivity becomes j v oc v^ 1 ^ and 
j u oc ^ _1 / 4 , respectively. 



3. Acceleration timescales 

The efficiency of shear acceleration can be studied 
by determining the acceleration timescales for typical 
shear flow velocity profiles. In the present context we 
will concentrate on well-collimated (cylindrical) rel- 
ativistic AGN jets, leaving a detailed analysis of ex- 
panding relativistic jet flows to the near future (Rieger 
& Duffy, in preparation) . In general, at least three dif- 
ferent shear scenarios may be distinguished by means 



of observational and theoretical arguments 



by n 



3.1. Gradual shear flow along the jet axis 

As noted in the introduction, there is growing ev- 
idence today for an internal velocity structure paral- 
lel to the jet axis, with the simplest scenario consist- 
ing of a fast moving inner spine and a slower moving 
boundary layer. In addition, Laing et al. (1999) have 
shown that the intensity and polarization systemat- 
ics in kpc-scale FR I jets are suggestive of a radially 
(continuously) decreasing velocity profile v z (r) . In or- 
der to estimate the particle acceleration efficiency in a 



longitudinal gradual shear flow, the following applica- 
tions seem thus to be particularly interesting (see also 
Fig.fJ): (i) a shear flow with velocity profile decreasing 
linearly from relativistic speed w Zimax to it Ziin i n over a 
characteristic length scale Ar = (f2 — ri), and (ii) a 
Gaussian- type velocity profile with core radius r c . 

3.1 .1 . Linearly decreasing velocity profile 

It can be shown (see Rieger & Duffy 2004) that for 
a linearly decreasing velocity profile the shear accel- 
eration timescale becomes 



^acc (0 



(Ar) 



A lb (r) 2 [l + lb (r) 2 u z (r) 2 /c 2 ] {Au z f 



(G) 



where A (7') is the mean free path for a particle with 



commoving Lorentz factor 7', Au z = (u z>r 







and 7fc(r) = 1/(1 — u z (r) 2 / c 2 ) 1 / 2 is the local bulk 
flow Lorentz factor. For nonrelativistic it Zim i n and 
7&(ri) larger than a few, the minimum acceleration 
timescales is thus of order 



iacc(n) 



3(Ar) 2 
7b(ri) 4 Ac 



(7) 



Note, however, that if M Zjm ; n is still relativistic with 
M z ,min ~ 0.7 u z , max as suggested, for example, in 
the case of FR I jets 0], the minimum accelera- 
tion timescale might easily be an order of magnitude 
larger. The evolution of the acceleration timescale 
as a function of r is shown in Fig. [21 The accelera- 
tion timescale increases significantly (up to a factor 
~ 2 • 10 3 for the application in Fig. |2K) while going 
outwards with r due to the decrease in flow velocity, 
suggesting that the higher energy emission will be lo- 
cated closer to n. 

Any particle acceleration process will have to com- 
pete at least with radiative energy losses. If one re- 
quires efficient shear acceleration, this thus translates 
into an upper limit for the allowed width Ar of the 
transition layer. In the case of synchrotron radia- 
tion, for example, the cooling timescale is given by 
^cooi oc to 3 7 ,_1 B~ 2 , with m the rest mass of the par- 
ticle. For a gyro-dependent particle mean free path, 
i.e., A ~ r g{l'), the condition iacc^i) < ^cooi thus 
implies an upper limit 



Ar £ (c 4 / e 5 / 2 ) 7fc (r 1 ) 2 m 2 B- 3 / 2 . 



(8) 



For typical pc-scale parameters, i.e., magnetic field 
strength of B ~ 0.01 Gauss and 76 (ri) ~ 10, one 
thus requires a width Ar e < 0.003 pc and Ar p < 
(m p /m e ) 2 Ar e (improving the estimate in |2l|) for ef- 
ficient electron and proton acceleration, respectively. 
While efficient acceleration of protons is thus very 
likely, efficient acceleration of electrons on the pc-scale 
and below appears quite restricted, i.e., only possible 
if the transition layer is relatively thin. Note, however, 
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Figure 1: Schematic illustration of a linearly decreasing and a Gaussian- type velocity profile for a shear flow parallel to 
the jet axis, with « z , ma x up to an inner radius r\ and u z , m in at r 2 , and r c — r\ for the Gaussian profile. 




Figure 2: The acceleration timescales for a gradual, parallel shear flow with linear decreasing velocity profile as a 
function of the radial coordinate. Left: Shear flow with bulk Lorentz factors 7f,(ri) = 5 and 7, respectively, r 2 = 50 ri 
and u z (r2) = 0.05it z (ri). The minimum acceleration timescale is defined by t acc (min) = i a cc(ri) Right: Shear flow with 
conditions more appropriate for FR I jets, i.e., bulk Lorentz factors 7f,(>"i) = 2 and 4, respectively, r 2 — 2r 1 and 
u z {r 2 ) = 0.7u*(n). 



that for powerful large-scale relativistic jets (such as 
the one in 3C 273) with typical parameters 7&(j"i) ~ 5 
and B ~ f0~ 5 G, only Ar e < 30 pc would be re- 
quired, suggesting that efficient electron acceleration 
may be more likelier in large-scale relativistic jets. 
Note also, that for A scaling with the gyro-radius, 
^acc °c icooi oc I/7', so that radiation losses are no 
longer able to stop the acceleration process once it 
has started to operate efficiently. 

3.1 .2. Gaussian velocity profile 

Consider for comparison a Gaussian-type velocity 
profile, which for r > 77 decreases as 

u z (r) = it z , min + (Au,) exp [-(r - r 1 f/{2 r 2 c )] , (9) 

i.e., asymptotically approaches u Zim - ln , where r c is the 
core radius and Au z the velocity difference as defined 



above. The acceleration timescale then becomes, cf. 
Rieger & Duffy (2004), 

t / n = t exp[(r - n) 2 /rg] 

(r/r c - rx/r c ) 2 j b {r) 2 [1 + 7 h (r) 2 u z (r) 2 /c 2 } 

(10) 

with to — 3cr^/[A (Au 2 ) 2 ] and r > r\. The minimum 
of the acceleration timescale is then no longer located 
at 77 , but still close-by. The evolution of the accelera- 
tion timescale as a function of the radial coordinate r 
is shown on the left hand side of Fig. [3] The acceler- 
ation timescale now first decreases up to a minimum 
before eventually increasing again significantly due to 
the diminishing velocity shear. An example, compar- 
ing the acceleration timescale for a Gaussian with the 
one for a linear-decreasing velocity profile, is shown 
on the right hand side of Fig. [21 illustrating the range 
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Figure 3: Left: The acceleration timescales for a longitudinal gradual shear flow with a Gaussian- type velocity profile 
as a function of the radial coordinate for "fb(ri) = 5, and tt z ,mm = 0.5 it z , max and tt z , m in = 0.05u Zimax , respectively, 
r c = n and r < 4 n . Right: The ratio of the acceleration timescale t^ cc for a Gaussian profile versus the acceleration 
timescale v£. c for a linear decreasing profile, shown using 7f,(ri) = 5, r c — n, r2 = 4ri, and w Zim in = 0.5u z , ma x- 



over which a Gaussian is slightly preferable. Note, 
however, that for the chosen parameters the minimum 
Gaussian acceleration timescale is still higher by a fac- 
tor of a few than the minimum acceleration timescale 
for the linear decreasing profile. 

3.2. Non-gradual shear flow along the jet 
axis 

Suppose a particle becomes so energetic that its 
mean free path A exceeds the width Ar of the velocity 
transition layer. The particle may then be regarded as 
essentially experiencing a non-gradual velocity shear, 
i.e., a discontinuous jump in the flow velocity. It has 
been pointed out by Ostrowski [Tfl[l7| that the jet side 
boundary in powerful AGN jets may naturally repre- 
sent a relativistic example where efficient non-gradual 
shear acceleration may occur. The minimum accel- 
eration timescale (measured in the ambient medium 
rest frame) can be estimated using Monte-Carlo sim- 
ulations 0, and appears to be comparable to the 
minimum timescale of nonrelativistic shock accelera- 
tion in the Bohm limit (A ~ r g ), i.e., 

r 

t acc ~ 10 — provided r g > Ar , (11) 

where r g denotes the particle gyro-radius. Balanc- 
ing the acceleration timescale with the synchrotron 
cooling timescale then translates into an upper limit 
for the maximum attainable particle Lorentz factor of 
7 < 4 • 10 7 (m/m e ) B~ x l 2 . Efficient acceleration thus 
requires 

Ar<2-1(T 8 (—) B^ 2 pc. (12) 
\m e J 

Hence, for a pc-scale magnetic field strength of B ~ 
0.01 G, a width of Ar e < 2 • 10~ 5 pc and Ar p < 



(m p /m e ) 2 Ar e would be needed in order to allow for 
electron and proton acceleration, respectively. Effi- 
cient electron acceleration on the pc-scale thus ap- 
pears unlikely, while efficient acceleration of protons 
could be well possible as long as the particle mean 
free path remains smaller than the width of the jet. 
Note that for the large-scale case with B ~ 10 -5 G 
one would still need Ar e ~ 0.6 pc, suggesting that 
electrons are also not accelerated efficiently on kpc- 
scales. 

3.3. Gradual transversal shear flow 

While the formation of relativistic jets generally 
seems to be connected to the presence of an accretion 
disk, a central black hole is not necessarily required 
If jets indeed originate as collimated disk winds, 
however, a significant amount of rotational energy of 
the disk may be channeled into the jet leading to a flow 
velocity field which is characterized by an additional 
rotational component (cf. |2fj , l2l| for more details). In 
order to evaluate the acceleration potential associated 
with such a shear flow, we have recently analyzed the 
transport of energetic particles in rotating relativistic 
jets 20] . Based on an analytical (mixed- frame) ap- 
proach, and starting with the relativistic Boltzmann 
equation and a (BKG) relaxation scattering term, we 
have studied the acceleration of particles in a cylin- 
drical jet model with relativistic outflow velocity u z 
and different azimuthal rotation profiles. It turned 
out, for example, that in the case of a simple (non- 
relativistic) azimuthal Keplerian profile, where shear 
effects dominate over centrifugal effects, local power 
law distributions f(p) oc p~( 3+Q ) are generated for 
t cx p a , a > 0, whereas for more complex rotational 
profiles (e.g., flat rotation, where centrifugal effects 
become relevant) steeper spectra may be produced. 
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In order to gain insights into the acceleration effi- 
ciency, we may consider a flow velocity field with con- 
stant relativistic u z and azimuthal Keplerian rotation 
profile of the form O(r) = fife (nn/f ) 3 ^ 2 , where f2/- is a 
constant. The acceleration timescale is then of order 
(see H3) 

' acc(r) * { lb (r)m-ui/ci)nl (i) ' (13) 

assuming the flow to be radially confined to a re- 
gion r- m < r < rj, where rj is the jet radius and 

7i,(r) = 1/[1 - ^(r) 2 r 2 /c 2 - ul/c 2 ] 1 ' 2 is the local flow 
Lorentz factor. According to Eq. (|13|1 efficient particle 
acceleration generally requires a region with signifi- 
cant rotation. As the acceleration timescale increases 
considerably with r, the higher energy emission will 
generally be concentrated closer to the axis (i.e., to- 
wards smaller radii). A comparison of acceleration 
and cooling timescales indicates that efficient electron 
acceleration on the pc-scale is only possible close to 
fin, whereas proton acceleration is not subject to such 
a restriction. Note again that for A oc 7, the ratio 
of acceleration timescale to cooling timescale is inde- 
pendent of 7, suggesting that losses are no longer able 
to stop the acceleration process once it has started to 
operate efficiently. 



• Gradual and non-gradual shear acceleration pro- 
cesses usually work very well for protons, which 
indicates that it may be possible to accelerate 
cosmic-rays to ultra-high e nerg ies along power- 
ful relativistic jets, cf. also |17lll8| . 

• Shear acceleration may generally represent an 
efficient way of dissipating a considerable part 
of the kinetic energy of the jet, thus possibly ac- 
counting for the substantial deceleration of jets 
in lower-power radio sources (e.g., 26]). 

• If a jet and its shear layer decelerates and ex- 
pands significantly with distance as expected, 
for example, in FR I sources, the minimum shear 
acceleration timescale will increase as well, sug- 
gesting the existence of a distance scale above 
which shear acceleration of electrons may no 
longer overcome the energy losses. 

• Compared with a simple homogeneous jet 
model, the emission properties of astrophysical 
shear flows may be much more complex, e.g., 
the inverse Compton emission can be strongly 
boosted as each jet part will see an enhanced 
radiation field from the other parts, cf. 



4. Conclusions 

There is growing evidence today that astrophysical 
jets may indeed be characterized by a significant shear 
within their flow velocity field. Internal jet rotation, 
for example, is likely to be present at least in the initial 
parts of the flow. On the other hand, a significant 
and continual velocity shear parallel to the jet axis, 
is expected for most powerful jet sources. Applying 
the results derived above thus suggests the following 
conclusions: 

• Gradual shear acceleration of electron occur- 
ring in small-scale relativistic jets may natu- 
rally account for a steady second population 
of synchrotron-emitting particles, contributing 
to the observed emission in addition to shock- 
accelerated [T^j ones. Depending on the level of 
turbulence (i.e., the momentum index a for the 
mean scattering time) shear-accelerated elec- 
trons may lead to very flat synchrotron emission 
spectra. On the other hand, if the Alfven veloc- 
ity is very high, the radio and optical emission 
properties of large-scale relativistic jets, such as 
the one in 3C 273 0, 0] , may be dominated by 
synchrotron emission from relativistic electrons 
accelerated via second-order Fermi processes (cf. 
|24|). The concrete details, however, will depend 
on the energy scale at which shear acceleration 
takes over |21|. 
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